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ABSTRACT 

We present the detection of a giant radio halo (GRH) in the Sunyaev-Zerdovich (SZ)- 
selected merging galaxy cluster ACT-CL J0256.5+0006 (z = 0.363), observed with 
the Giant Metrewave Radio Telescope at 325 MHz and 610 MHz. We find this cluster 
to host a faint (Se io = 5.6 ± 1.4 mJy) radio halo with an angular extent of 2.6 arcmin, 
corresponding to 0.8 Mpc at the cluster redshift, qualifying it as a GRH. J0256 is one 
of the lowest-mass systems, M 5 qo,sz = (5.0 =b 1 . 2 ) x 10 14 M 0 , found to host a GRH. We 
measure the GRH at lower significance at 325 MHz (S 325 = 10.3±5.3 mJy), obtaining 
a spectral index measurement of = 1.0to. 9 - This result is consistent with the 
mean spectral index of the population of typical radio halos, a = 1.2 ± 0 . 2 . Adopting 
the latter value, we determine a 1.4 GHz radio power of Pi. 4 gh z = (1-0 ± 0.3) x 10 24 
W Hz -1 , placing this cluster within the scatter of known scaling relations. Various 
lines of evidence, including the ICM morphology, suggest that ACT-CL J0256.5±0006 
is composed of two subclusters. We determine a merger mass ratio of 7:4, and a line- 
of-sight velocity difference of v± = 1880 ± 210 km s _1 . We construct a simple merger 
model to infer relevant time-scales in the merger. From its location on the Pi. 4 gh z -Lx 
scaling relation, we infer that we observe ACT-CL J0256.5±0006 just before first core 
crossing. 

Key words: Galaxies: clusters: individual (ACT-CL J0256.5±0006) - Galaxies: clus¬ 
ters: intracluster medium - radio continuum 
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1 INTRODUCTION 


Multiwavelength observations of galaxy clusters provide a 
wealth of information about the physics of the intraclus¬ 
ter medium (ICM) and its relationship with cluster galax¬ 
ies. The optical and X-ray bands have historically been 
used to identify merger activity via optical substructure 


fCarter &; Metcalfe 

1980; 

Geller &; Beers |1982| Rhee & 

Katgert 1987; Dress 

er & S 

lectman: 1988; Rhee, van Haar- 


lem fe Katgert| |1991| |Wen &; Han 2013| ) and morpho- 
logical parameters determined from X-ray images (|Mohr, 
Fabricant &; Gellerj |1993| JJeltema et ah] |2005| |Q’Hara 
et al. | 2006 1 jSantos et ah |2008| ). In the last decade, a 
link has been found between a cluster’s merger status 
and the presence of large-scale diffuse synchrotron emis¬ 
sion (see |Brunetti fe~J ones 2014, and references therein). 
This cluster-scale radio emission, dubbed a giant radio 
halo (GRH) if ~Mpc in size, exhibits a steep spectrum 
and has no obvious link to the individual cluster galax¬ 
ies citepBuote.2001.GRH, FerettiGiovannini.2008.GRH, 
Ferrari.2008.GRHReview, Feretti.2012.Review. Radio halos 
(RHs) appear to trace the non-thermal ICM and typically 
have spectral indices of a ~ 1.1-1.5. However, ultra-steep 
spectrum radio halos (USSRHs, a ~ 1.6-1.9), presumably 
associated with more pronounced synchrotron ageing, have 
also been detected within the population (|Brunetti et al. 


2008 Dallacasa et al.|2009 Venturi et al.|20LL 

The existence of USSRHs is predicted by one of the cur¬ 


rent leading theories for the origin of RHs (Brunetti et al. 


2008), namely the turbulent re-acceleration model in which 
the synchrotron emission is powered by turbulence gener¬ 
ated during cluster mergers ( Brunetti et al.|200l Petrosian 


2001 Brunetti V Lazarian||2011 Beresnyak et al. 2013). In 

this model one expects an USSRH to be seen when the tur¬ 
bulent energy in the cluster has decreased sufficiently for it 
to be less efficient in accelerating high energy electrons in the 
cluster. This scenario can also explain the observed bimodal¬ 
ity in scaling relations between the 1.4 GHz RH power and 
thermal cluster properties, in which clusters are observed 
to be either radio loud or radio quiet. This dichotomy has 
been observed in cluster samples selected via X-ray lumi¬ 
nosity ( Brunetti et al.|[2007 Cassano et al. 2008) and the 
Sunyaev-Zel’dovich (SZ) effect ( Sunyaev V Zel’dovich|1972 ), 
although it is less pronounced in the latter case (Sommer 
V Basu||2014 ). In practice, one anticipates a population of 
clusters in transition between these two states that will have 
intermediate radio power. 

The observed bimodality was initially thought to be due 
in part to selection effects in the cluster sample ( Basu|20l2 ), 


but this has since been ruled out (Sommer & Basu 2014 
Cuciti et al. 2015). A more likely reason is a physical effect 


related to the cluster evolutionary state. Magnetohydrody¬ 


namic (MHD) simulations by Donnert et al. (2013) show 


that a RH is a transient phenomenon that exhibits a rise 
and fall in radio halo emission over the course of a merger. 
This evolutionary model suggests that for a merging clus¬ 
ter, the observable diffuse radio emission depends strongly 
on the phase of the merger in which the cluster is being 
observed, which likely contributes to the scatter in the ob¬ 


served Pi. 4 GHz scaling relations with thermal cluster prop¬ 
erties. 

Moreover, one would expect to find two separate types 
of systems that populate the intermediate region of radio 
power: late-stage mergers with old RHs that are in the pro¬ 
cess of switching off, and early-stage mergers in which the 
radio halo emission has recently switched on but not yet 
reached its maximum radio power. The former scenario is 
a possible explanation for some of the observed USSRHs, 
which are starting to fill in the region between the correla¬ 
tion and upper limits. Clusters that are in the early stages 
of merging would also be interesting systems to identify and 
study as they would complete the evolutionary picture; how¬ 
ever due to their expected low radio power, they are poten¬ 
tially more difficult to detect. 


In line with the above, Cassano et al. (2010) find that 


the observed dichotomy is strongly related to cluster dy¬ 
namical state, with morphologically disturbed systems host¬ 
ing RHs. However, several RH non-detections in merging 
clusters are seemingly incongruent with this trend (A141, 
A2631, MACSJ 2228: |Cassano et al.|2010IA119 : |Giovannini| 
V FerettT||2000 and A21 46: |Russell et al.|20lf ). In the case 
of A2146, Russell et al. (2011) postulate that the lack of a 


RH in this strongly-merging system is due to the relatively 
low mass of the cluster. They estimate a radio power upper 
limit more than an order of magnitude below the correlation. 
Low-mass systems are expected to generate less turbulent 
energy during their mergers, yielding weaker synchrotron 
emission, and hence RHs that are too faint to observe with 


current telescopes. The era of LOFAR (Vermeulen 


SKA precursors such as MeerKAT (Booth & Jonas 


2012 ) 


2012 ) 


and ASKAP ( DeBoer et al.|2009 ), and the SKA itself (Tay- 


lor||2013 ) will bring with it highly sensitive observations of 


these systems, and should reveal the underlying RH emis- 


In this paper we present the detection of a GRH in a 
low-mass system that we argue is in the early stages of merg¬ 
ing. As discussed, such early-stage merging systems are in¬ 
teresting because they allow us to probe the full evolutionary 
cycle of GRHs and are expected to fill in the intermediate 
region in radio halo power. 

The paper is organised as follows: we present existing 
multi wavelength data on ACT-CL J0256.5+0006 in Section 
[2] and we describe the radio observations and data reduc¬ 
tion process in Section [3] with the radio results presented in 
Section [4] X-ray and optical morphological analyses are dis¬ 
cussed in Sections |5.1| and |5.2| respectively. We construct a 
model for the merger geometry in Section [6] and infer merger 
time-scales from this model in Section |71 We conclude with 
a discussion in Section [8] In this paper we adopt a ACDM 
flat cosmology with Ho — 70 km s -1 Mpc -1 , — 0.27 and 

Ua s= 0.73. In this cosmology, at the redshift of our cluster 
(z=0.363), one arcminute corresponds to 305.8 kpc. We as¬ 
sume S u oc throughout the paper, where S u is the flux 
density at frequency v and a is the spectral index. Colour 
versions of all figures are available in the online journal. 


2 ACT-CL J0256.5+0006 

ACT-CL J0256.5+0006 (hereafter J0256) lies at z=0.363 
and was detected by the Atacama Cosmology Telescope 
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Table 1. Published properties of J0256. 


R.A. (hh mm ss.s) 

02 56 33.0 a 

Dec. (dd mm ss.s) 

+00 06 26.3 a 

redshift 

0.363 b 

L 5 oo,x (10 44 ergs s -1 ) 

3.01 ± 0.36 c 

U 500 (10 4 arcmin 2 ) 

3.4 ± 1.0 d 

M 500 ,x (10 14 M 0 ) 

c 9+1-1 e 
°- Z -0.9 

M 50 o,sz (10 14 M 0 ) 

5.0 ± 1.2 d 


R.A. and Dec. (J2000) of the SZ peak of the cluster, with an 
astrometric accuracy of 5-10+ 


Menanteau et al. 


(2013) 


c Integrated 0.1-2.4 keV X-ray luminosity using the 
spectrum of |Majerowicz et al.||2004| , corrected for the 
cosmology adopted in this paper. 

d Integrated Compton y-parameter and B12 SZ mass from 
|Hasselh eld et al. |2013| . 

e Total mass for the main cluster component using /3-model fit 
“a” for the NE region ( |Majerowicz et al.[2004| . 


(ACT; Kosowsky 2006) equatorial SZ cluster survey with a 
148 GHz decrement signal-to-noise ratio of 5.4 (Hasselfield 
et al.|2013 ). It was first identified in ROSAT PSPC data and 
is included in the Bright SHARC catalog (RX J0256.5+0006; 
Burke et al.||1997 ). Majerowicz et al. ( 2004| identify J0256 

as undergoing a major merger based on observations carried 
out with XMM-Newton. 

In the following sub-sections we describe the exist¬ 
ing multiwavelength data for J0256 in the X-ray ( XMM- 
Newton ), optical (Gemini), millimetre (ACT), and radio 
(VLA) bands. The relevant cluster properties are given in 
Tabled] 


2.1 X-ray 


Majerowicz et al. (2004), hereafter M04, carry out a compre¬ 


hensive X-ray study of J0256 based on their 25.3 ks XMM- 
Newton observations (obs ID: 00560230 lQ The X-ray im¬ 
age shows two components in the direction of the cluster: a 
bright main component and a less luminous structure to the 
west. To investigate whether these are physically connected 
or serendipitously aligned, M04 fit an elliptical /3-model to 
the hot gas distribution of the main component, excluding 
point sources and the western component. After subtraction 
of the best-fit model from the data, the residuals reveal that 
the western component is a small galaxy cluster exhibiting a 
comet-like morphology, with the tail to the west (see Figure 

2 in M04). This orientation indicates that gas in the sub¬ 
cluster is undergoing ram pressure stripping as it interacts 
with the main cluster component. Based on the orientation 
of the subcluster isophots away from the main component 


and numerical simulations by Ricker X Sarazin (2001), M04 


conclude that the subcluster has not yet passed through the 
main cluster centre and thus that J0256 is in the pre-core 
crossing stage of its merger. 

For the full cluster, M04 determine a temperature of 


1 The XMM-Newton observations include European Photon 
Imaging Camera (EPIC) data from the two MOS (Metal Oxide 
Semi-conductor) CCD arrays and the pn CCD array. 


T — 4 . 9 lS^ keV within ~ O.8R500 and a bolometric X-ray 
luminosit^Jj of Lx = (7.88 ± 0.53) x 10 44 erg s -1 , which 
is over-luminous compared to the Ly^-T relation measured 
by Arnaud X Evrard (1999). M04 conclude that this dis¬ 


crepancy between observed and predicted luminosity, cou¬ 
pled with their evidence for ram pressure stripping of the 
subcluster, suggests J0256 is not in dynamical equilibrium. 
Using XSPEc[^]to model the M04 spectrum using the cosmol¬ 
ogy in this paper, we determined a 0.1-2.4 keV band-limited 
luminosity of L 5 oo,x = (3.01 ±0.36) x 10 44 erg s -1 , incorpo¬ 
rating a conservative 10 % uncertainty due to the spectrum 
being extracted within ~ 0 . 8 i? 5 oo- 

From /3-model fitting, M04 calculate an X-ray mass for 
the main cluster component of M 500 ,x = 3.71q;| x 10 14 M q , 
which is equivalent to M 5 oo,x = 5 . 2^9 x 10 14 M o using 
the cosmology in this paper. The M04 M 200 mass for the 
main cluster component is M 200 ~ 5.9 X 1O 14 M 0 . Using 
count rates in the residual map in the region of the subclus¬ 
ter and translating this into a luminosity, they estimate the 
M 200 mass of the subcluster to be between 1-2 xlO 14 M 0 
and determine a merger mass ratio of ~ 3:1. However, this 
calculation requires several broad assumptions due to a lack 
of ancillary data, making the result somewhat uncertain. 


2.2 Millimetre 

Wide area, untargeted SZ surveys detect large numbers of 
galaxy clusters via inverse Compton scattering of cosmic mi¬ 
crowave background (CMB) photons by electrons within the 
hot ICM, which causes a distortion of the CMB spectrum 
in the direction of clusters. ACT is a 6 m telescope that ob¬ 


serves the millimetre sky with arcminute resolution (Swetz 
et al.|[ 20 TT ). Between 2008 and 2011 , ACT surveyed a 455 
deg * 2 * strip centred at 5 = -55°, as well as a 504 deg 2 strip 
centred at 5 — 0° overlapping the Sloan Digital Sky Survey 


(SDSS) Stripe 82 (Marriage et al. 2011 Hasselfield et al. 
2013), at 148, 218, and 277 GHz. ACT has detected over 


ninety clusters via the SZ effect. 

J0256 was identified in the ACT equatorial 148 GHz 
map, with a decrement signal-to-noise ratio of 5.4 for a fil¬ 
ter scale of $500 = 7.06' (see Hasselfield et al. 2013 here¬ 
after H13, for details). H13 investigated prescriptions for 
the pressure profile used to obtain a Y500—M500 scaling 
relation, where Y+o is the integrated Compton parame¬ 
ter. HI3 investigated several profiles computed from sim¬ 
ulations (e.g., Battaglia et al. 2012) or empirical models 
(e.g., Arnaud et al.|2010 ), leading to a SZ mass range of 2.9 


xlO 14 M 0 < M500 < 7.5 x 10 14 M q for J0256, taking into ac¬ 
count the range of uncertainties on all mass estimates. The 
pressure profile from |Battaglia et"ah (2012) is currently pre¬ 
ferred, and in this paper we use the corresponding SZ mass 
estimate of M 5 oo,sz — (5.0 ± 1.2) x 1O 14 M 0 . 


2.3 Optical 

The ACT collaboration has completed spectroscopic obser¬ 
vations of J0256 using Gemini and identified 78 cluster mem¬ 
bers (Sifon et al. 2015). This distribution of spectroscopi- 


2 Corrected for the cosmology used in this paper. 

3 https://heasarc.gsfc.nasa.gov/xanadu/xspec/ 
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Figure 1 . Cluster member galaxies with spectroscopic redshifts from Gemini identified on an SDSS r-band image. Blue boxes (red 
circles) denote members with higher (lower) redshifts than the systemic cluster redshift of z = 0.363. Large, bold symbols mark the 
BCGs of both kinematic components. The 148 GHz Compton y SZ contours are superposed. The contours start at a level of 2.0 X 10 -5 , 
increasing towards the centre in steps of 1.25 X 10 -5 . The black X marks the cluster SZ peak. 


cally confirmed cluster members is ~85% complete within 
R 200 , up to an r-band magnitude limit of 22. Using this red- 
shift information, we can estimate an independent dynami¬ 
cal mass and re-examine the merger geometry proposed by 
M04 (see Section |5^2| below) . The cluster members are shown 
in Figure [l] where red circles (blue boxes) denote members 
that are at lower (higher) redshifts than the cluster redshift 
of z — 0.363. We identify these two sets of galaxies as sep¬ 
arate kinematic components (see Section [ 5 T 2 ] below), each 
of which has a brightest cluster galaxy (BCG) that is indi¬ 
cated by a large, bold symbol. If the cluster is not in the 
core passage phase of its merger, the superposition of the 
two populations in the plane of the sky indicates that the 
merger is occurring at least partially along the line-of-sight. 


2.4 Radio 


J0256 has been mapped at 1.4 GHz in the NRAO VLA Sky 
Survey (NVSS; [Condon et al. 1998) and the Faint Images 
of the Radio Sky at Twenty-Centimetres (FIRST; |Becker, 


White &; Helfand||1995 ) survey and at 74 MHz in the VLA 


Low-Frequency Sky Survey (VLSS; Cohen et al.|2007 ). Fig¬ 
ure [2] shows the cluster region in each of the three sky sur¬ 
veys. Only one point source is detected in the 1.4 GHz 
survey data at R.A. and Dec. (J2000) of 02h56m34s and 
+00d065m03. Its NVSS and FIRST fluxes are 4.8 =L 0.4 


mJy and 3.66 d= 0.27 mJy, respectively. This source is not 
detected in the VLSS data; however, there is a source 1.16 / 
away, closer to the SZ peak of the cluster, detected 3cr above 
the map noise. The rms and resolution of each image is given 
in the caption for Figure [2] 


3 NEW RADIO OBSERVATIONS 

We observed J0256 with the Giant Metrewave Radio Tele¬ 
scope (GMRT) as part of an ongoing project involving the 
radio follow-up of ACT equatorial clusters. Initial observa¬ 
tions were carried out for 10 hours at 610 MHz in August 
2012 (PI: Knowles), using a 33 MHz bandwidth split into 
256 channels and a 16s integration time. The data were ac¬ 
quired in the polarization channels RR and LL, and the 
total on-source time was 7.5 hrs. Flux and bandpass cal¬ 
ibrator 3C48 was observed at the beginning, middle, and 
end of the observing block. This source was also used to es¬ 
timate the instrument’s antenna gains and ionospheric phase 
calibration which in turn were used to correct observations 
of the target field. A second set of 8-hour observations was 
carried out at 325 MHz on the GMRT using Director’s Dis¬ 
cretionary Time (PI: Knowles) in July 2014. This dataset 
has a central frequency of 323 MHz with a total bandwidth 
of 33 MHz made up of 256 channels and an integration time 
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Figure 2. Postage stamp images of the J0256 cluster region at 1.4 GHz from NVSS (left) and FIRST (middle), and at 74 MHz from 
VLSS (right). The dashed black circle denotes R 500 centred on the SZ peak, which is marked by a black X. The image resolutions, from 
left to right, are 40" X 40", 6.4" x 5.4", and 75" x75". The rms is given in the upper right corner and the beam is indicated by the 
yellow ellipse at lower left in each image. The colour scales are all in units of mJy beam -1 . 


Table 2. GMRT observations. 


Frequency® 

Observing date 

On-source 

Integration 

Bandwidth^ 

^synth? P- a -° 

rms noise c 

HPBW 

n d 

"max 

(MHz) 


time (hrs) 

time (s) 

(MHz) 

("x", °) 

(pjy beam -1 ) 

(arcmin) 

(arcmin) 

610 

Aug 2012 

7.5 

16.1 

29.1 

5.7 x 4.3, 71.3 

26 

43 ± 3 

17 

325 

Jul 2014 

6.5 

8.1 

31.2 

9.8 x 8.2, 76.1 

72 

81 ± 4 

32 


a Observing frequency. 
b Bandwidth remaining after flagging. 

c Synthesised beam and rms noise of the full-resolution images, where p.a. denotes the beam position angle. 
d Maximum recovered scale. 


of 8s. The total on-source time was 6.5 hrs. As with the 610 
MHz observations, 3C48 was used as the sole calibrator. Ob¬ 
servational details are given in Table [2] The pointing centre 
for both sets of observations was the same and was defined 
to be that of the SZ peak, given in Table [l] 

The 610 MHz and 325 MHz data were subjected to the 
same calibration procedure, which is based on AIPS (NRAO 
Astronomical Image Processing System), SPAM ( |Intema 
et al.|[2009| ), and Obit ( Cotton|[2008| tools. The main cali¬ 
bration steps are outlined here. First, strong radio frequency 
interference (RFI) is removed by statistical outlier flagging 
tools. As a compromise between imaging speed and spectral 
resolution losses due to bandwidth smearing, the datasets 
are then averaged down to 24 channels. Phase calibration 


starts from a model derived from the VLSS (Cohen et al. 


2007) and the NVSS (Condon et al. 1998), followed by a 


succession of self-calibration loops. To compensate for the 
non-coplanarity of the array, we use the polyhedron (facet- 
based) wide-field imaging technique available in AIPS. We 
perform several rounds of imaging and self-calibration, in¬ 
specting the residual visibilities for more accurate removal of 
low-level RFI using Obit. To correct for ionospheric effects, 
we then apply SPAM calibration and imaging. The presence 
of strong sources in the field of view enables one to derive 
direction-dependent (DD) gains for each source and to use 
these gains to fit a time variable phase screen over the entire 


array. The phase screen was used during imaging to correct 
the full field of view for ionospheric phase effects. 

As J0256 lies at close to zero declination, bright sources 
in the field are subject to strong north-south sidelobes that 
interfere with emission in the cluster region. To reduce the 
impact of these bright sources during further imaging, we 
modeled and subtracted all sources in the field outside of 
a 13 arcminute radius centred on the cluster, leaving a 
dataset with only the inner portion of the field. This edited 
n^-dataset was then imported into the Common Astron¬ 
omy Software Applications package (CASA; [McM ullin et al. 


2007) for imaging. 


For each dataset we created several target field images 
all with Briggs robust R = 0 weighting (Briggs] 1995). We 
first made full resolution (FR) images, shown in Figures Al 
(610 MHz) and |A3| (325 MHz) in the Appendix, using all of 
the iru-data, cleaning until the residuals were noise-like. We 
then created high-resolution (HR) images in the following 
way. As the 610 MHz data have more long baselines than 
the 325 MHz data, we matched the r^-coverage of the two 
datasets by selecting a iru-range from 4 kA (~52") to 30 
kA (~6 X/ ), and imaging using a 25 kA outer taper. The HR 
images were cleaned until their residuals showed no indica¬ 
tion of emission in the cluster region. The clean components 
from the HR images were used as compact source models 
and were subtracted from the m>-data to create a point 
source-subtracted datasets. Using these datasets, we imaged 
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Table 3. Properties of the different radio images created. Values 
in brackets are for the 325 MHz images when different from the 
corresponding 610 MHz images. 


Image ID 

^min 

(aremin) 

0 * 
Umax 

(aremin) 

Point sources 

removed 

FR 

0.08 (0.13) 

17 (32) 

No 

HR 

0.13 

0.86 

No 

PSSUB-FR 

0.08 (0.13) 

17 (32) 

Yes 

PSSUB-LR 

0.84 

17 (32) 

Yes 

LR* 

1.30 (1.26) 

17 (32) 

Yes 


* The highest resolution available, defined by the synthesised 
beam. 

* The largest scale to which the image is sensitive, defined by the 
shortest baseline/wn-wavelength. 

* PSSUB-LR convolved with a V Gaussian. V corresponds to 
-3.5 kA. 



36.0 2:56:32.0 36.0 2:56:32.0 


Right ascension Right ascension 

0 12 3 4 

Figure 3. Left : 610 MHz full-resolution (FR) image of the clus¬ 
ter region showing seven radio sources. Right : 610 MHz full 
resolution image of the same region after subtracting the point 
source model from the uu-data (PSSUB-FR). The colour scale is 
in units of mjy beam -1 and is the same for both panels. 


at full resolution (PSSUB-FR) to visually check that the 
point source subtraction was successful. 610 MHz HR and 
PSSUB-FR images of the cluster region are compared in the 
left and right panels of Figure [3] respectively. The PSSUB- 
FR image shows no visual indication of residual emission 
from the compact sources; however, we nevertheless investi¬ 
gate contamination from the source removal process in Sec¬ 
tion [472] Once satisfied, we re-imaged with a uv -cut of < 4 
kA and an outer taper of 3 kA to gain sensitivity to diffuse 
emission on scales of 1 Mpc, creating point source sub¬ 
tracted, low-resolution (PSSUB-LR) images. We convolved 
each PSSUB-LR image with a 1' Gaussian, providing better 
sensitivity to extended features while retaining useful data, 
to create our final smoothed, point source subtracted, low- 
resolution (LR) maps shown in Figures A2 (610 MHz) and 


A4 (325 MHz) in the Appendix. The final LR 610 MHz (325 


MHz) map has a maximum angular resolution of 17' (32 7 ) . 
A summary of the different images created is given in Table 

0 


4 RADIO RESULTS 

With the angular resolution and short baselines of the 
GMRT, we are able to investigate emission from both com¬ 
pact sources and extended diffuse structures. In the follow¬ 
ing, we discuss our results from both the 610 MHz and the 
325 MHz datasets. 


4.1 Compact radio sources 

There are seven bright radio sources in the cluster region 
identified in both 325 MHz and 610 MHz full-resolution 
maps, five of which are associated with spectroscopically 
confirmed cluster members. The 610 MHz HR contours can 
be seen in the left panel of Figure [4] along with source labels. 
The only source detected in NVSS and FIRST, as discussed 
in Section [2M] is detected in our maps as S7. The flux densi¬ 
ties and spectral index we measure for this source, provided 
in Table imply a consistent 1.4 GHz flux density of 4.61 
=b 0.64 mjy. 

Several of these sources exhibit resolved tail emission, 
possibly due to merging activity in the cluster. The BCG of 


the subcluster is associated with the radio source S5. This 
source has a wide extension to the west of the galaxy, and 
although our highest resolution image cannot resolve finer 
structure within the extended tail, it may be a bent narrow 
angle tail radio galaxy contorted by ram pressure stripping 
due to the merger (Bl iton et al.|1998 ). The multi-frequency 
radio properties of all seven sources are given in Table [4] 
Here and in Section |4.3.2| the spectral indices are deter¬ 
mined using a Monte-Carlo simulation, in which we draw 
from Gaussian flux density distributions with means and 
widths represented by the flux densities and their uncer¬ 
tainties, respectively. The spectral index and uncertainties 
are then determined from the median and 68th percentiles 
of the resulting spectral index distribution. 


4.2 Point source contamination 

To unveil any low surface brightness extended cluster emis¬ 
sion, the HR radio sources, particularly in the cluster region, 
have to be removed from the 'uu-data as described in Section 
[3] Although the point source removal is reasonably success¬ 
ful, as is clear from the right panel of Figure [3] it is not exact. 
In order to quantify the residual (low) level of contamina¬ 
tion, we perform a statistical analysis of the LR image using 
both radio source and random off-source positions in the 
following way: 

(i) In the HR image, we select a large number (>100) of 
random off-source positions. 

(ii) For each position, we calculate the LR map flux den¬ 
sity in a LR beam-sized area centred on that position. 

(iii) From this set of flux densities we calculate the mean, 
/Rand, and standard deviation, cr ran d, of the distribution. We 
expect /irand to be close to zero for Gaussian noise. 

(iv) We then select all sources outside of the cluster re¬ 
gion that are detected above 5 a in the HR map; we find 28 
resolved and 53 unresolved sources. 

(v) We repeat steps (ii)-(iii), now using the point source 
positions, /i ptS rcs quantifies the bias in subtraction of point 
source emission. cr ptS rcs contains both the map uncertainty 

and a measure of the noise added by the subtraction process, 

2 _ 2,2 

CRyst? 1-6- CTp tsrcs ^rand ' ^syst • 
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Table 4. Properties of cluster region radio sources. Source labels are shown in the left panel of Figure The given R.A. and Dec. 
are for the peak source emission in the 610 MHz map. Flux errors include 10% measurement uncertainties. The uncertainties on a are 
determined via numerical methods, as described in Section [tT] 


Source 

RA 

(hms) 

DEC 

(dms) 

Type 0 

S^io 

(mJy) 

8325 

(mJy) 

a» 

Notes 

SI 

02 56 35.5 

00 06 11.0 

c 

0.56 ± 0.08 

0.69 ± 0.12 

0.33 ± 0.31 


S2 

02 56 35.9 

00 06 27.9 

T 

2.17 ± 0.24 

3.32 ± 0.37 

0.67 ± 0.21 


S3 

02 56 33.8 

00 06 28.8 

C 

2.17 ± 0.24 

3.76 ± 0.41 

0.87 ± 0.21 

associated with BCG of main component 

S4 

02 56 32.6 

00 06 30.9 

T 

1.20 ± 0.15 

1.93 ± 0.23 

0.75 ± 0.24 


S5 

02 56 30.4 

00 06 01.8 

T 

4.14 ± 0.43 

9.71 ± 0.98 

1.35 ± 0.19 

associated with BCG of subcluster 

S6 

02 56 32.2 

00 05 50.8 

C 

0.42 ± 0.08 

0.59 ± 0.12 

0.54 ± 0.38 

foreground source 

S7 

02 56 33.8 

00 05 02.0 

T 

7.71 ± 0.78 

11.39 ± 1.15 

0.62 ± 0.20 

detected in NYSS and FIRST C 


a C: compact; T: resolved with tailed emission. 

b Spectral index between 325 MHz and 610 MHz (S u oc v ~ a ). Errors are determined via Monte Carlo methods (see text for details). 
c Extrapolating S6io to 1.4 GHz using a S7 gives S 1400 — 4.61 d= 0.64 mJy, which is consistent with the values quoted in Section [2.4| 
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Figure 4. Left : GMRT 610 MHz high-resolution (6.5" X 5.0 // , p.a. 78.9°) contours of the J0256 emission, overlaid on the SDSS gri-band 
image. The high-resolution (HR) image lcr noise level is 31 pjy beam -1 and the contours are [3,10,20,40,80] x lcr. The HR beam is shown 
as the yellow ellipse in the lower left corner. Individual radio galaxies are labelled from SI to S7. Flux densities for these sources can be 
found in Table The red X marks the position of the SZ peak. Right : Smoothed XMM-Newton MOS X-ray contours (arbitrary levels 
from the reprocessed image starting at 12 counts per second per square degree - see Figure [8] in section |5T| ) , overlaid on the smoothed 
low-resolution (LR) 610 MHz image of the GRH in J0256. The LR radio image is obtained after subtracting the compact source emission 
from sources SI to S7 (positions marked by black crosses). The red circle ( blue square) indicates the BCG of the main (subcluster) 
component. The positions of the BCGs coincide with the X-ray peaks of each component. The LR synthesised beam (79.6 // X 76.8", p.a. 
-86.9°) is shown as the yellow ellipse. The purple X marks the position of the SZ peak. The radio colour scale has units of Jy beam -1 . 


The results of this analysis are given in Table [5] We 
find that we are systematically over-subtracting a low level 
of point source emission, more so when the sources are re¬ 
solved. Moreover, the subtraction process does add a small 
but non-negligible amount of noise into the LR image, as ex¬ 
pected. Using the relation in step (v) above, this systematic 
noise is cr sys t,6io = 0.3 mJy beam“^ in the 610 MHz map 
and (j syst ,325 = 1.0 mJy beam“^ in the 325 MHz map. We 
incorporate these systematic and random residuals into our 
final flux density measurements (see Section 4.3.1). 


A graphical representation of this process is shown in 
Figure [5] In the HR and LR maps, we stack on the source 
and random off-source positions separately. The left panels 


of Figure [5] show the stacked results from the HR map. As 
expected, the random positions produce a noise-like result 
and the stacked source positions produce a clear compact 
source at the centre. 

Repeating this process in the LR image, we find a neg¬ 
ative stacked signal slightly off-centre from the source posi¬ 
tion, in agreement with the over-subtraction implied by in 
Table [5] The shifted peak is due to the varying noise in the 
map, shown by the random stacked result (middle panels 
of Figure [5|. We note that the rms of the LR source and 
off-source stacked maps are comparable. 

As a final check, we stack on the radio source posi¬ 
tions in the PSSUB-FR map and smooth this result to the 
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Table 5. Results of the systematic and statistical tests to quan¬ 
tify the residual point source contamination in the low-resolution 
maps. All values are in units of mJy beam“^. 


V 

(MHz) 

Quantity 

Source Positions 
Compact Resolved 

All 

Random 

Positions 

Number of sources 

53 

28 

81 

116 

610 


-0.075 

-0.082 

-0.077 

0.013 


(7 

0.547 

0.822 

0.655 

0.586 

325 

n 

-1.073 

-1.920 

-0.971 

0.273 


a 

3.109 

2.470 

2.693 

2.503 


—0.077 ± 0.073 mJy beam”^, leading to a 5% larger cor¬ 
rected flux density for the halo. However at 325 MHz, 
M 325 , P tsrcs = —0.971 d= 0.299 mJy beanWj, which is a bias 
measured at a significance of 3cr that leads to a fractional 
flux density increase of over 50%. We thus correct the mea¬ 
sured flux densities and incorporate the systematic uncer¬ 
tainties introduced by the point source removal into the flux 
density uncertainties. We also include ~10% absolute flux 
calibration and residual amplitude errors 
. The final flux density, S u 
Aare calculated as follows: 

Su — *S7,meas (AG,ptsrcs X Ns) (1) 

A Si = (0.1 Su) 2 + (<tLb + <st) x (N S ) (2) 


Bhatnagarj 2004 
ing uncertainty, 


Chandra, Ray & 
, and correspond- 


Table 6. GRH properties. Subscripts denote frequencies in MHz 
unless otherwise stated. 


S6io (mJy) 5.6 ± 1.4 

5 325 (mJy) 10.3 ± 5.3 

<a 610 1 n+ 0 - 7 

^325 -L.U q g 

Pi.4GHz (10 24 W Hz-^t l.o ± 0.3 
LLS6io (Mpc)* * 0.8 


1 Extrapolated from 5610 using a spectral index of a = 1.2 d= 0.2. 

* Largest linear size of the GRH, corresponding to 2.6'. 


where cr rms is the central map noise, cr sys t is the systematic 
error due to point source removal, and Ns is the number 
of independent beams within the flux aperture. We measure 
integrated halo flux densities of $610 = 3.6 ±1.4 mJy and 
5*325 = 10.3 ± 5.3 mJy. The additional contributions to the 
flux density uncertainty lower the significance of the 610 
MHz detection to 4cr which is low, but still reliable. The 325 
MHz flux, however, now has a signal-to-noise of less than 2. 
Higher sensitivity observations at 325 MHz are required to 
reliably confirm our detection at this frequency. 


same resolution as the LR map. These results are shown 
in the right panels of Figure [5] There is a net residual af¬ 
ter source subtraction mostly caused by imperfect subtrac¬ 
tion of resolved sources, the peak of which is ~ 10 % of the 
peak brightness of the average source in the stacked HR 
map. When we smooth to the same beam as the LR map 
(lower, right panel), we largely recover the structure of the 
LR stacked source result (upper, middle panel). 


4.3 Diffuse emission 

After removal of the radio sources in the field, the LR 610 
MHz map, shown in Figure |A2| reveals distinct extended 
emission in the cluster region with a 6a peak above the map 
noise. The 3 a angular extent of the emission is 2.6 / , corre¬ 
sponding to a physical scale and largest linear size (LLS) 
of 0.8 Mpc at the cluster redshift. Due to the centralised 
position and size of this emission, we classify it as a giant 
radio halo, making J0256 one of the lowest-mass clusters to 
host one known to date. The right panel of Figure [4] shows 
the 610 MHz GRH overlaid with smoothed X-ray contours. 
The GRH roughly follows the X-ray emission and is centred 
on the cluster SZ peak. The GRH radio properties are listed 
in Table [5] Our LR 325 MHz map is shown in Figure |A4| 
The radio peak of the GRH lies to the west of the cluster 
SZ peak and is only marginally detected at a significance of 
3cr in the 325 MHz map. 


4-3.1 Flux measurements 

The flux density is measured within an aperture of ra¬ 
dius 90 7/ , centred on the 610 MHz emission such that all 
610 MHz halo flux is captured. From the results of the 
point source contamination analysis in Section [4. 2 | the bias 
at 610 MHz is only at the la level, i.e., /i6io,ptsrcs = 


4-3.2 Spectral index 

We can estimate a theoretical spectral index for the GRH in 
J0256 from the distribution of measured radio halo spectral 
indices from the literature, shown in Figure [ 6 ] Assuming this 
cluster is in the early stages of merging, based on the X-ray 
morphology determined by M04 (see Section [2 .1 1 above), we 
expect J0256 to host a young, and therefore flatter spectrum, 
radio halo. We therefore exclude the USSRHs (a ^ 1.6) from 
the literature and use the mean and rms of the remaining 
17 radio halo spectral indices to determine our theoretical 
value and error respectively. We determine a spectral index 
for the typical radio halo population of a = 1.2 ± 0 . 2 . 

Our measured spectral index, 25 = l-Oj^g, obtained 
using 56io and the noisy S325 measurement, is consistent 
with the above value. However, given the large uncertainties 
on (* 325 , driven by the large error on S325 , we choose to adopt 
the spectral index of the regular radio halo population, a = 
1.2 ± 0.2, to extrapolate our measured GRH flux density to 
other frequencies. 


4-3.3 Radio power 


The 1.4 GHz GRH radio power, Pi. 4 GHz is correlated with 


thermal cluster properties and cluster mass (Cassano et al. 
2013). To constrain Pi. 4 GHz, we use our 610 MHz flux den¬ 


sity measurement and the assumed spectral index from the 
previous section to extrapolate a flux density at 1.4 GHz. We 
account for the effect of redshift on the flux density and ap¬ 
ply a k-correction to calculate a halo radio power of Pi. 4 GHz 
at (1.0 ± 0.3) xlO 24 W Hz -1 in the cluster rest frame. The 
error on Pi. 4 GHz is propagated from the spectral index un¬ 
certainties. We note that the radio power is consistent with 
the non-detections in NVSS, FIRST, and VLSS, as it cor¬ 
responds to a GRH surface brightness far below the noise 
levels of these surveys. 
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arcseconds arcseconds arcseconds 


Figure 5. Results from stacking on radio source and random off-source positions in the 610 MHz maps, where all colour scales are in 
units of mjy beam -1 and the map axes are in arcseconds. All maps are centred on source positions. Left panels: Stacked images in the 
HR map using radio source (upper) and off-source (lower) positions. The elliptical beam is 6" X 5". Middle panels: Same as the left 
panels, but for the smoothed, source-subtracted, low-resolution (LR) map. The beam here is 80" X 70". Right panels: Radio source 
stacked maps from the PSSUB-FR image (upper) smoothed to the LR beam (lower). 



Figure 6. Distribution of all measured GRH spectral indices 
in the literature. The bulk of the values are taken from IFeretti 
|et al.|(|2012| > with three updated measurements from |Venturi et aT 
(|2013| > and new GRHs from |Bonafede et al.| ( |2014a| ) and Bonafedej 
|et al.| |2014b| ). USSRHs (a ^ 1.6) are shown in light grey. 


J0256 is shown as the red star on the radio power cor¬ 
relations in Figure [7] The cluster lies within the scatter, and 
on the same side, of all three correlations from the literature. 
J0256 appears to lie slightly further away from the Piaguz~ 
Y500 relation, compared to its relative position in the other 
planes. However, the position of the cluster, relative to the 
distance away from each correlation, is consistent within the 
error bars for the cluster mass and thermal parameters. 


5 CLUSTER MORPHOLOGY 

As current observations favour a theory of merger-driven 
radio halo formation, it is important to understand the dy¬ 
namical state of J0256. With the X-ray and optical redshift 
information available to us, we can perform a morphological 
analysis of J0256. 


5.1 X-ray morphology 


Visual inspection of the reprocessed J0256 X-ray image in 
Figure[8]shows the cluster to be disturbed. This image is pro¬ 
duced by following the ESAS reduction thread for extended 
X-ray source^] in which the A1 and Si lines are modelled 
in XSPEC. The image has been both exposure-corrected and 
background-subtracted. We exclude the pn camera data as 
the pn CCD is marginally flared and has a chip gap near the 
cluster core. Point sources were removed during the repro¬ 
cessing using the cheese task in the standard XMM SAS 
tookQ The masked regions are shown in Figure [8]by yellow, 
dashed, excluded circles. 

In order to quantify the level of substructure in the re¬ 
processed XMM-Newton combined MOS1 + MOS2 image, 


we follow the work of Cassano et al. (2010) by calculating 


4 http://heasarc.gsfc.nasa.gov/docs/xmm/esas/esasimage/ 
esasimage_thread-str.html 

° http://xmmssc-www.star.le.ac.uk/SAS/xmmsas_20121219_ 
1645/doc/cheese/index.html 
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Figure 8. Combined 0.5-2.0 keV MOS 1 +MOS 2 XMM-Newton 
image of ACT-CL J0256.5+0006 with the 610 MHz radio halo 
contours overlaid. The contours start at 3cr = 0.001 Jy beam -1 , 
and increase in steps of la. The X-ray image has been exposure- 
corrected and background-subtracted. The region masked after 
removing point sources is indicated by the yellow, dashed, ex¬ 
cluded circle. The image is unbinned and has been smoothed by 
a Gaussian with a kernel radius of 6 pixels (1 px = 2.5 // ). The 
units of the colour scale are counts per second per square degree. 


three morphological parameters. To determine the measure¬ 
ment uncertainty on each of our parameters, we adopt the 
simulation method of Bohringer et al. (2010) whereby a Pois¬ 
son resampled X-ray image is used to compute the standard 
deviation of a parameter measurement, which is then used 
to estimate the measurement uncertainty. 


5.1.1 Concentration parameter, cs b 


The concentration parameter, proposed by |Santos et al.| 
(2008) as a probe of cluster substructure, is the ratio of the 
cluster core and the larger-scale X-ray surface brightnesses. 
We calculate the concentration parameter as 


S(< 100 kpc) 
S(< 500 kpc) ’ 


(3) 


where S is the X-ray surface brightness within a particular 
radius, centred on the X-ray peak. We determine a value of 
c SB = 0.151 =b 0.007 for J0256. 


5.1.2 Centroid shift, w 


Figure 7. Radio halo detections and upper limits from the lit¬ 
erature showing correlations between the 1.4 GHz radio power 
and cluster thermal parameters — (a) Pi. 4 vs Lx, (b) Pi. 4 vs 
l+oo 5 and (c) Pi. 4 vs M 500 , SZ- Black solid (open) circles and grey 
open triangles are giant radio halos (USSRHs) and upper lim¬ 
its, respectively, from Cassano et al. (2013), with recent GRHs in 
PLCK147.3-16.6 ( |van Weeren et al.||2014| and El Gordo ( |Lind-| 
|ner et ah][2 014) shown as blue squares. The position of J0256 is 
shown as a red star. The best fit to the GRH detections and asso¬ 
ciated 95% confidence interval is from |Cassano et al. (2013]) and 
are shown by the black line and grey shaded region, respectively. 


Poole et al. (2006) show that, compared to other X-ray mor¬ 


phological estimators, the centroid shift is the most sensitive 
to cluster dynamical state and least sensitive to cluster im¬ 
age noise. It is defined as the rms deviation of the projected 
separation between the X-ray peak and the centre of mass 
in units of the aperture radius, R ap , computed in a series 
of concentric circular apertures centred on the cluster X-ray 


peak (Mohr, Fabricant &; Geller 


- - . -1I9931.-- . 

Vaughan et al.||2008 l. Following |Cassano et al. (2010), the 

iperture radius is decreased in steps of 5% from a maximum 


O’Hara et al. 2006 
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aperture of radius R ap = 500 kpc to 0.05 R ap . We compute 
the centroid shift as 


N 


^ T ^(A i -(A }) 2 


1 1/2 


X /? ’ ^ 
-I+a P 


where A * is the distance between the X-ray peak and the 
centroid of the zth aperture. We measure a value of w — 
0.045 ± 0.006 for J0256. 


5.1.3 Power ratio, P 3 /P 0 


The power ratio of a cluster is calculated using a multipole 
decomposition of the potential of the two-dimensional pro¬ 
jected mass distribution. The idea of using the power ratio 
of the X-ray surface brightness to probe the underlying mass 


distribution was first introduced by Buote & Tsai (1995) and 


has since been widely used as an indication of substructure 


within a cluster (Jeltema et al.|2005 

Ventimiglia et al. 2008 

Bohringer et al.|2010 Cassano et al. 

2010). We use the nor- 


malised hexapole moment, P 3 /P 0 , which is the lowest power 
ratio moment providing a clear measure of substructure 
(Bohringer et ah] 2010), calculated in an aperture of radius 
R ap = 500 kpc centred on the X-ray cluster centroid. For 
J0256, we calculate a value of P 3 /P 0 = (3.76 ±0.30) x 10 -6 . 


Redshift 

0.350 0.355 0.360 0.365 0.370 0.375 



Velocity (km/s) 


Figure 9. Histogram showing the redshift distribution for 78 
spectroscopically confirmed cluster members. Here v = 0 is de¬ 
fined as the cluster systemic redshift of z = 0.363, and the bin 
width is 420 km s —1 . A bimodal fit of two Gaussians (dashed 
black) is shown along with the constituent main component (thick 
red; /a = 0.361 ± 0.001,(7 = 0.004 ± 0.001) and sub-component 
(thin blue; g = 0.369 ± 0.002, a = 0.003 ± 0.001). The verti¬ 
cal thick red (thin blue) dashed line shows the velocity of the 
BCG for the main (subcluster) component. A single Gaussian fit 
(/ 1 = 0.363 ±0.002, <7 = 0.005 ±0.001) is shown by the dot-dashed 
black curve. 


5.1.4 Comparison with the literature 


Using the methods described in Sections |5.1.1|5.1.3| jCas-| 


sano et al. (2010) study the morphological parameters for 


all clusters in the GMRT Radio Halo Survey (Venturi et al. 


2007 2008) and find a link between cluster dynamical state 


and the presence of a radio halo. They define a cluster to 
be dynamically disturbed if its morphological parameters 
satisfy the following conditions: csb < 0.2, w > 0.012 
and P3/P0 > 1.2 X 10 -7 . The majority of dynamically 
disturbed clusters are found to show radio halo emission. 
All of the parameter values we determine in our analysis 
of J0256 ( csb = 0.151 ± 0.007, w = 0.045 ± 0.006 and 
P3/P0 = (3.76 ± 0.30) x 10 -6 ) satisfy the above conditions 
for a merging cluster, as expected. 


We note that the Cassano et al. (2010) results were ob¬ 


tained using Chandra data whereas our results are obtained 
with XMM-Newton data, which has a larger PSF. To in¬ 
vestigate the effect of the different instruments on the vari¬ 
ous morphological parameters, we use archival Chandra and 
XMM-Newton data on a known merging cluster, A2631, and 
compare the derived morphological parameters from each 
image. The exposure times for the Chandra and XMM ob¬ 
servations, after flare rejections, are 16.8 ks and 13.4 ks, 
respectively. We find that the resolution difference between 
the two instruments has a negligible effect on the concentra¬ 
tion or centroid shift parameters. However, the power ratio 
is higher in the XMM image by a factor of between 2 and 
5, depending on the level of smoothing applied. Even with 
a reduction by a factor of five, the P 3 /P 0 value for J0256 
(7.5 xlO -6 ) is still well above the threshold of 1.2 x 10 -7 
for dynamically disturbed clusters. 


5.2 Optical redshift distribution 

X-ray morphological parameters are largely insensitive to 
substructure along the line of sight. To gauge any disturbed 
morphology in this direction, we use the redshift distribution 
of 78 spectroscopically confirmed cluster member galaxies 
(see Section [T3| above). This distribution is shown in Fig¬ 
ure [5] there is an indication of bimodal structure in the 
histogram. 


5.2.1 Statistical analysis using GMM 


To gauge its significance, we perform a Gaussian mixture 
model (GMM) analysis of the member galaxy redshifts. We 
use the GMM code developed by Muratov X Gnedin (2010) 


to fit a 2-mode Gaussian mixture to our data and compare 
it to a unimodal fit. The code calculates the kurtosis of the 
distribution, K , and the maximum log likelihood, logL, to 
which each model converges. For a bimodal fit, the peak 
separation of the modes relative to their widths, D , is also 
calculated. A statistically significant bimodality would have 
K < 0, D > 2, and a log-likelihood value greater than that 
for a unimodal fit. Parametric bootstrapping of the uni¬ 
modal distribution is performed to determine the probabili¬ 
ties of the observed K , D , and logL difference values being 
sampled from a unimodal distribution. The latter probabil¬ 
ity defines the confidence interval at which a unimodal fit 
can be rejected. 

The results of our analysis are given in Table [7] The 
multi-variance bimodal mixture model and unimodal Gaus¬ 
sian fits are superimposed on the distribution in Figure [9] 
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Table 7. GMM statistics from the redshift distribution of 78 cluster members. All errors are at the lcr level. 


kurtosis, K 
peak separation, D 

-0.260 

2.64 ± 0.82 






Distribution type 

n 

Statistics 

y cr 2 

logit 

Bootstrapping (%)t 

K D x 2 

Unimodal 

78 

0.363 ± 0.001 

0.005 ± 0.000 

299.6 

- 

- 

Bimodal, multi-variance 

53.9 ± 15.9 
24.1 ± 15.9 

0.360 ± 0.002 
0.369 ± 0.003 

0.004 ± 0.001 

0.003 ± 0.001 

300.7 

49.0 

46.6 69.4 


t The maximum log likelihood to which the fit converges. The difference in log L values defines a x 2 proxy. 
$ Measure of how likely it is that the same statistic can be drawn from a unimodal model. 


shown by the dashed and dot-dashed curves respectively. 
The data satisfy the K < 0 and D > 2 criteria for bi¬ 
modality, with the largest log L value coming from the multi¬ 
variance bimodal fit. The improvement in the logL value for 
the multi-variance bimodal model relative to the unimodal 
model is not significant due to the difference in degrees of 
freedom; a likelihood-ratio test indicates that the bimodal 
fit is rejected in favour of the unimodal fit at 53%. According 
to the parametric bootstrapping, the unimodal distribution 
is consistent with the data at the 69% level when only the 
log(L) probability is considered, with bootstrapped proba¬ 
bilities of K — 49% for the kurtosis, and D — 47% for the 
peak separation. A unimodal fit thus cannot be ruled out. 

However, statistical tests run on mock bimodal 
datasets, with the same population ratio and number of 
members as our real data, show that the GMM test re¬ 
sults in a log(L) probability of 70% or higher about 10% of 
the time. An unambiguous bimodal preference is only con¬ 
sistently achieved once the total population size is greater 
than 200. This implies that, when the distribution size is 
small, the GMM test could show a slight preference for a 
unimodal fit even when the input redshift distribution is bi¬ 
modal, given the population ratio of our true sample. There¬ 
fore, with the available number of redshifts for J0256, the 
GMM log(L) test is not a strong discriminator between the 
two models. However, based on the following additional ev¬ 
idence, we adopt the bimodal model in further analysis of 
this cluster. 


Firstly, there are two BCGs (cluster members with the 
lowest SDSS magnitudes) that are spatially separated, as 
seen from the SDSS image in Figure [l] which are also sepa¬ 
rated in velocity space as shown in Figure [9] This provides 
support for the existence of two distinct galaxy populations. 
These galaxies coincide with the peaks in the XMM-Newton 
X-ray emission (see Figure [4]). Secondly, the DS test, which 
measures the deviation of the velocity distribution in spa¬ 
tially localised regions of a cluster with respect to the cluster 
as a whole, indicates the presence of substructure in J0256, 


with Sa < 0.01 at the 68% confidence level (Sifon et al. 


2015). 


We use the GMM code to provide, for each member 
galaxy, the probability that the galaxy belongs to each of the 
kinematic components in the multi-variance bimodal case. In 
the following section we use these probabilities to calculate 
physical properties for the cluster and its components. 


5.2.2 Velocity dispersions and dynamical masses 


By fitting a 2-mode GMM to our data, each cluster member 
is assigned a probability of belonging to each of the modes. 
These probabilities can be used to determine the mean and 
variance for each mode by integrating over ah members and 
weighting by the probabilities. Since we have a discrete num¬ 
ber of member galaxies, the mean and variance for compo¬ 
nent n are given by 


Z n = ( Z) n 


'l2iPn( z i) 


(5) 


= (( z -z) 2 ) = ^ Pn(Zl)Z|2 -/ z f 

U j J2 iP n(zi) {ZK ' 


( 6 ) 


where n £ {1, 2}, z% is the redshift of the i-th member galaxy, 
and Pn(zi) is the probability that this member belongs to the 
n-th component. The mean and variance of each mode in the 
redshift distribution correspond to the peak redshift and ve¬ 
locity dispersion for each kinematic component, respectively. 
We use the velocity dispersion and the galaxies-based scal¬ 


ing relation from Munari et al. (2013) to determine M 200 


and i7 2 oo for each componenFl using a value of h m 0.7 in 


the Munari et al. (2013) equation. Using the concentration 


parameter from Duffy et al. (2008), we integrate a NFW 


profile (Navarro, Frenk &; White 1997) and interpolate to 
determine M 500 and R $00 • The results are given in Table [8] 
with all uncertainties determined via bootstrapping. We fol¬ 
low the same process using the unimodal fit, the difference 
being that the probability for every member is 1. 

From the mean redshifts of the components, we find 
a line-of-sight velocity difference of v± = 1880 ±210 km 
s _1 . We also calculate individual component masses of 
M 500 ,main — (3.23 ± 0.66) x 10 14 M® and M 500 ,subcl. — 
(1.83 ± 0.74) x 10 14 M q , leading to a merger mass ratio of 
7:4, smaller than but within the errors of the ~ 3:1 ratio 
determined by M04. Combining the component masses, we 
calculate a cluster dynamical mass of M 5 oo,o P t = (5.06 ± 
0.99) x 10 14 M©, which agrees with the SZ cluster mass given 
in Table [l] to better than 0.5cr. The combined M 2 oo,opt mass, 
M 2 oo,o P t = (7.66 ± 1.54) x 10 14 M Q , agrees within lcr with 
the estimated M 2 oo,x total cluster mass range from M04 of 
M 2 oo,x ~ 9.7— 11.1 x 10 14 M©, assuming a 15% uncertainty 
on their M 2 oo,ne valu<Q 


6 M 200 = (47r/3)p200^1oO 

7 Corrected for the cosmology in this paper. 
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Table 8. Optical statistics of the two cluster components from 78 spectroscopic galaxy redshifts. v pec is relative to z — 0.363. 


Component No. of galaxies 

^mean 

'Cpec 

<7 

AF 200 

M 500 

#200 

#500 



(km s -1 ) 

(km s -1 ) 

(1O 14 M 0 ) 

(10 14 M 0 ) 

(Mpc) 

(Mpc) 

main cluster 59 

0.361 ± 0.001 

-490 ± 100 

850 ± 70 

4.90 ± 1.03 

3.23 ± 0.66 

1.45 ± 0.11 

0.92 ± 0.06 

subcluster 19 

0.369 ± 0.002 

1390 ± 180 

690 ± 120 

2.76 ± 1.14 

1.83 ± 0.74 

1.20 ± 0.19 

0.76 ± 0.12 



Figure 10. Merger geometry of J0256 as per Majero wicz et ah] 
(2004). The small black dot represents the centre of the main 
cluster component and the red circle represents the centre of the 
subcluster. d m i n and d are the projected distance and physical 
distance between the two component centres, respectively. v± is 
the line-of-sight infall velocity and 0 is the impact angle. 


If we model the cluster as a single component, we esti¬ 
mate a total mass -ATsoOjtot — (7.74±0.02) x 10 14 M g , which 
is 2.3cr away from the SZ mass. This reinforces our argu¬ 
ment in favour of the bimodal model. The corresponding 
M 200 measurement, Af 2 oo,tot ~ 11.7 x 10 14 M©, still agrees 
with the total X-ray mass estimate from M04, although this 
comparison is not particularly meaningful given the large 
uncertainties on their estimate. 


6 MERGER GEOMETRY 


M04 construct a simple merger model for J0256 using pro¬ 
jected distances and the line-of-sight velocity difference be¬ 
tween the main and subcluster components. We adopt a sim¬ 
ilar approach but update two aspects: we use a more cur¬ 
rent cosmology and the increased number of galaxy spectro¬ 
scopic redshifts (78 vs. 4) discussed in Section 


2.3 


The op¬ 
tical galaxy redshift distribution also allows us to determine 
dynamical masses for the main and subcluster components. 

For simplicity, we assume the same merger geometry 
as in M04, schematically outlined in Figure [lO] Working 
in the rest frame of the main component, we assume the 
same simplification of a point mass subcluster and ignore 
dynamical friction. However, rather than using a /3-model, 
we assume the mass distribution of the main component is 
defined by a NFW profile (Navarro, Frenk & White][l997): 


M(< R) — AtvpqRI 


ln(l + c) - 


1 + c 


(7) 


where R s = R/c is a characteristic scale radius, c is the 
concentration parameter for radius #, and po is the typical 
NFW dark matter density for the cluster. Using the c(M, z) 
relation from Duffy et al. (2008) to determine c for our clus¬ 


ter, we have c = 3.018 and po = 5.497 x 10 M© Mpc 



v and cluster component separation d. The intersections of the 
two relations give the two possible solutions for v and d. The 
grey dotted lines and shaded regions indicate the solution for 
case one: v\ = 1930^^q km s -1 and d\ = FlOS^i^ 53 Mpc. The 
green dotted lines and shaded areas indicate the solution for case 
two: V 2 = 2640^gJ° km s -1 and d 2 = 0.338^q‘q2q Mpc. 



Figure 12. The line-of-sight velocity v± versus the component 
separation d (black, thick, solid) using the relation in equation [ 9 ] 
with the infall velocity v(d) given by equation [8] The solid thin 
red curve is the impact angle 6 as a function of d. The vertical and 
horizontal blue dashed lines indicate the values of d m i n = 237.6 
kpc and v± = 1880 ±210 km s -1 respectively, with the horizontal 
blue dot-dashed lines indicating the lower and upper limits for v ±. 
Figure E3 showed the two merger geometry solutions. Here the 
grey dotted lines and shaded areas indicate the d and 6 values and 
uncertainties for case one: d± = 1.11^!^ 35 Mpc and 6 \ = 12° 

The green dotted lines and shaded regions indicate the same, but 
for case two: cfo = 0.34 +q'q 2 Mpc and 62 = 45°^®. 
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Table 9. Merger geometry and time-scales from today for two possible cases with d m i n = 237.6 kpc and uj_ = 1880 ± 210 km s 1 . 



V 

d 

e 

-t A a 

tB b 

tc c 

r d 


(km s" 1 ) 

(kpc) 

(degrees) 

(Gyr) 

(Gyr) 

(Gyr) 

(%) 

case 1 

1930lJ®o 

1110+ 350 

iiiu_ 240 

12±| 

1 06 +0 ‘ 23 
l.UO — 0 20 

n 4e+ 0 - 21 
U -^°-0.22 

1 QQ + 0 - 20 

i.yy _0 22 

oc + 7 

oO — ig 

case 2 

26401“° 

340±2o 

45 1 ° 

1 41+ 0 - 03 
-*-•^-0.04 

0 12+ 0 - 03 
U - ±Z -0.03 

1 64 +0 ' 03 

± - Dq: -0.03 

46±2 


a Time since first virial crossing. 
b Time until core passage. 
c Time until second virial crossing. 

d Measure of how far along in the merger the cluster currently is, T = |tA/^tot| =1 |^a/(^C — £a)|- 


Using the above mass profile and modelling the grav¬ 
itational infall of the subcluster, we obtain the following 
relation between subcluster infall velocity, v, and physical 
separation, d, between the centres of the subcluster and the 
main component: 


ln(l + d/R s ) ln(l + c) 

d/Rs c J ’ 

( 8 ) 

where Mo = AirpoR^. The subcluster redshift z su b is greater 
than that of the main cluster component, z m am- As we ar¬ 
gued in Section H3 the X-ray emission pattern indicates 
that the subcluster is moving towards the main component. 
This implies that the impact angle must be less than 90°. 
Using simple trigonometry, it follows from the merger geom¬ 
etry in Figure [To] that 

OU+tv)*- 1 - < 9 > 


v\d) = 


2 GM 200 2GMq 


7^200 


R s 


where d m i n is the projected separation between the main 
component and the subcluster, and v± is the velocity differ¬ 
ence along the line-of-sight. 

Using the X-ray peaks of each component, the projected 
separation between cluster components is ~0.78 7 , which cor¬ 
responds to a physical projected distance of d m i n = 237.6 
kpc (as compared to 350 kpc in M04). In Section 5.2.1 we 
found v± = 1880 d= 210 km s -1 which is consistent with the 
value estimated by M04. Based on the X-ray arguments in 
Section |2.1| the two cluster components have begun inter¬ 
acting and we can place the following limits on the physical 
separation and the infall velocity: d m i n < d < R 200 and 
v > v±, where R 200 is the cluster radius for the main com¬ 
ponent. 

Simultaneously solving equations [8] and [9] with these 
constraints provides two sets of solutions for the merger 
model. These are listed in Table|9] with the graphical so¬ 
lutions given in Figures El an d |12| The uncertainties on v, 
d , and 9 are shown in Figures 11 (v and d) and [l2| (d and 
0), and are propagated from the uncertainties on the R 200 
mass and radius of the main cluster component, the R 200 
radius of the subcluster, and measured line-of-sight velocity 
difference. We consider these solutions in the next section 
to estimate relevant time-scales in the merger. 


7 MERGER AND RADIO HALO 
TIME-SCALES 

To better understand the formation history and mecha¬ 
nism (s) of GRHs, we would like to relate the GRH formation 
time-scale to the merger time-scale. It is possible to model 
the physics of turbulent re-acceleration using simulations. 
Donnert et al. ( 2013| ) (hereafter D13) used MHD simulations 
of a 10 lt> M© and 8:1 merger to study the strength and pat¬ 
tern of diffuse radio emission at various merger stages. They 
found that the cluster needs to have been actively merging 
for a minimum amount of time, approximately 15% into the 
merger, such that there is sufficient turbulence generated, 
before the radio emission switches on. 


7.1 Estimates for merger time-scales 


To estimate the merger time-scales for J0256 we assume a 
simple merger taking place in a linear fashion along the 
merger axis determined by the impact angle, 0 , schemati¬ 
cally outlined in Figure [13] In Section |2.1[ we ruled out a 
scenario in which the subcluster has already passed through 
the core. In Figure pffi] we isolate three distinct times during 
the merger: (A) first virial crossing; (B) core passage; and 
(C) second virial crossing. Even though we refer to virial 
crossing, we use R 200 as a proxy for the viria l radius. 

From the optical analysis in Section 


5.2 


Mpc and Rloo cl — 1-20 Mpc. First virial crossing thus oc¬ 
curs when the centres of the two components are initially 
2.65 Mpc apart. The distances associated with the three 
merger stages are Ra — 2.65 Mpc — d, Rb = d , and 
Rc = 2.65 Mpc + d , where d is the current physical sep¬ 
aration for the two model solutions listed in Table [9] 

For each merger model solution found in the previous 
section, we compute the merger times 

r R n dR ' 


tn — 


rttr 

Jro 


r(R') 


( 10 ) 


where n G {A, B, C}, Ro is the observed position of the sub¬ 
cluster, and v NFW is the velocity function given in equation 
[8] The total time of the merger, at least for the first passage, 
is given by Uot — tc —t a- We define the relative time phase 
of the merger as the ratio Y m |tA/Uot|- The results for each 
model solution are given in Table [9] 

For case 2, we find that J0256 would have completed 
first virial crossing 1.41 Gyr ago with ~120 Myr until first 
core passage occurs. This puts the cluster T 2 = 46t 2 % °f 
the way into its merger. In case 1, J0256 is closer to the be¬ 
ginning of its merger with ~460 Myr until first core passage. 
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Figure 13. Schematic showing the relative position of the subcluster (red circles) to the main cluster (black circle) overlaid on the X-ray 
image at three different times during the merger: (A) first virial crossing; (B) core passage; and (C) second virial crossing. The centre of 
the main cluster is marked with a white cross while the centre of the subcluster at each interval of the merger is shown by a red diamond. 
The blue cross and dashed circle denotes the current position of the subcluster. The dashed black line represents the merger axis and 
d m in is the projected distance between the two cluster components. All circles denote R 200 of the respective components. Even though 
we refer to virial crossing, we use R 200 as a proxy for the virial radius. 


The time-scales for case 1 result in J0256 having a relative 
time phase of Ti = 35^ 8 %. According to D13, these con¬ 
clusions lead to very different theoretical predictions for the 
observed strength and morphology of the radio emission. In 
the following section we compare our time-scale results with 
the D13 simulations. 

7.2 Comparison with MHD simulations 

The simulated radio powers and morphologies in D13 are 
for observations at 1.4 GHz of a massive 1O 15 M 0 cluster 
undergoing a 8:1 mass ratio, plane-of-the-sky merger. J0256 
is about 50% of the total simulated mass but has a much 
smaller mass ratio of 7:4. As the strength, and hence observ¬ 
ability, of the radio emission is related to cluster mass and 
the amount of turbulent energy created during a merger, 
we caution that, for the specific case of J0256, the following 
comparison with the D13 results can at best be qualitative 
due to the above differences between J0256 and the sim¬ 
ulated cluster. MHD simulations for the particular case of 
J0256 would be required for a more accurate comparison. 

To compare our merger time-scales with the MHD sim¬ 
ulations of D13 we need to convert our values into their 
time frame. From the X-ray snapshots of their simulated 
merger (see their Figure 3), we estimate first and second 
virial crossings to occur at 0 Gyr and 2.56 Gyr respectively, 
giving £tot,Di3 = 2.56 Gyr, similar to the total merger time 
of 3.05 Gyr for J0256. Scaling our T values to this time-scale 
allows us to extrapolate expected radio power and general 
emission morphology for each case in Table [9] using the D13 
simulation. 

In case 1 we have Ti = 35^ 8 %, corresponding to 
£a,di3 = 0.90iQ 4g Gyr. Here, not enough turbulence is 
being generated to drive the diffuse radio emission and 
only compact radio source emission is observable in Fig¬ 
ure 3 of D13. Case 2 gives T 2 = 46^2%? corresponding to 
£a,di3 — 1.181q q 5 Gyr. Here the GRH is in the early stages 
of having switched on, according to Figure 3 of D13, and is 
gaining power. The X-ray image of J0256 shown in Figure [8] 


is a close visual match with the second panel of Figure 3 in 
D13, which has a relative time-scale similar to that of case 
2. This consistency is in contrast to case 1, where no diffuse 
radio emission is observable and the expected radio power 
lies in the realm of the upper limits on the Pi. 4 GHz-+x scal¬ 
ing relation. Thus our case 2 appears to be the more likely 
of the two merger geometry solutions for J0256: we observe 
what is likely a young radio halo. 


8 CONCLUSION 

We have detected a low surface brightness giant radio halo 
(-0.8 Mpc) in ACT-CL J0256.5+0006 with the GMRT at 
610 MHz, and obtained a marginal detection at 325 MHz. 
With an SZ mass of M 50 o = (5.0 ± 1.2) x 1O 14 M 0 , J0256 
is one of the lowest mass clusters currently known to host 
such emission. 

We measure halo flux densities of $610 = 5.6 + 1.4 mJy 
and S 325 = 10.3 + 5.3 mJy, giving a measured spectral index 
0:325 = l.O^o . 9 - Due to the unreliability of the 325 MHz 
measurements, we calculate a bandwidth- and k-corrected 
1.4 GHz radio power of Pi. 4 GHz = (1.0± 0.3) x 10 24 W Hz -1 
by extrapolating our 610 MHz flux density to 1.4 GHz using 
a theoretically motivated spectral index of a = 1.2 + 0.2. As 
the detection at 610 MHz is not highly significant, we do not 
draw strong conclusions about the radio morphology, but we 
do note that it roughly follows the thermal gas as seen in 
the X-rays and is centred on the cluster SZ peak. More data 
at 325 MHz would be required to confirm our detection at 
this frequency and obtain a more accurate measured spectral 
index. 

Using the X-ray and optical information available to us, 
we have investigated the morphology of J0256, concluding 
that this system consists of a main cluster component with 
an in-falling subcluster slightly in front and to the west of 
it. The merger mass ratio determined via new spectroscopic 
galaxy member redshifts is roughly 7:4, making it a major 
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merger event. We estimate a line-of-sight velocity difference 
between the two components of v± = 1880 d= 210 km s _1 . 

Using this information and assuming an NFW mass pro¬ 
file and a simple merger geometry defined by v, d, and 6 , we 
hnd two possible solutions for the merger time-scale. Defin¬ 
ing the merger time phase, T, to be the percentage of the first 
passage (between first and second virial crossings) already 
completed, we hnd that J0256 has a merger time phase of 
r 2 = 46^2% or Ti = 35ti 8 %. We compare these values with 
MHD simulations from Donnert et al. (2013) and conclude 
that J0256 is most likely ~47% of the way into its merger, 
with only ~100 Myr until first core passage. As the strength 
of the synchrotron emission is related to the amount of tur¬ 
bulent energy produced during a merger, a population of 
simulations varying in cluster mass and merger ratio would 
be useful in investigating the GRH formation rate for a wider 
range of models. 

Our discovery of a GRH in J0256 may help to provide 
some insight into whether GRHs exist in all merging clusters 
and whether the non-detections in known merging systems 
are due to a combination of a low-mass cluster and insuffi¬ 
cient sensitivity to diffuse emission, rather than to a com¬ 
plete lack of GRHs. More systems like J0256 will probe the 
full evolving population of GRHs, in particular the early- 
stage mergers, and potentially fill in the gap between radio 
upper limits and USSRHs in the Piagh.z—Lx plane, ft would 
be interesting to carry out a similar merger time-scale analy¬ 
sis for existing GRHs to probe the scatter in the radio power 
scaling relations. 
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APPENDIX A: FULL-RESOLUTION AND 
LOW-RESOLUTION RADIO MAPS 


In this appendix we provide the inner 30' x 30' of the full 
resolution and smoothed low resolution maps for both 610 
MHz and 325 MHz. In each image, the dashed circle ind i- 
cates the cluster scale $500 = 3.T from Hasse ffield et al.| 
(2013), centred on the SZ cluster peak, which is shown as a 
red or white X. The solid circle shows the 13 x radius outside 
of which we removed all compact emission before further 
imaging in CASA, as described in Section [3] 
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Figure Al. Inner 30' X 30' of the full-resolution (FR) 610 MHz map. The beam is 5.7" X 4.1" at p.a. 71.3°, and the map noise is a 
= 26 p,Jy beam -1 . The dashed black circle represents $500 = 3.1', centred on the cluster SZ peak shown by the red X. The 13' radius is 
shown by the solid black circle. 
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Right ascension 

Figure A2. Inner 30 r X 30' of the 610 MHz map. Greyscale is the low-resolution (LR), H-smoothed image. Red contours are the 
high-resolution (HR) [6, 20, 80]xlcr contours where la = 31 juJy beam -1 . The X and black solid and dashed circles are as in Figure | aT| 
The LR beam is 79.6 " X 76.8" at p.a. -86.9° and is shown by the blue ellipse in the lower left corner. The lcr noise in the LR greyscale 
image is 0.36 mJy beam -1 . 
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Figure A3. Inner 30' X 30' of the full-resolution (FR) 325 MHz map. The beam is 9.7" X 7.9" at p.a. 74.1° and the map noise is a 
77 /iJy beam -1 . The X and black solid and dashed circles are as in Figure [Al] 
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Figure A4. Inner 30' X 30' of the 325 MHz map. Greyscale is the low-resolution (LR), l'-smoothed image. Red contours are the 
high-resolution (HR) [6, 20, 80]xlcr contours where lcr = 71 /xJy beam -1 . The X and black solid and dashed circles are as in Figure [ aT| 
The LR beam is 79.4'' x 73.1" at p.a. 56.7° and is shown by the blue ellipse in the lower left corner. The lcr noise in the LR greyscale 
image is 1.18 mJy beam -1 . 













